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Abstract 

We link the study of positive quantum maps, block positive operators, 
and entanglement witnesses with problems related to multivariate polyno- 
mials. For instance, we show how indecomposable block positive operators 
relate to biquadratic forms that are not sums of squares. Although the 
general problem of describing the set of positive maps remains open, in 
some particular cases we solve the corresponding polynomial inequalities 
and obtain explicit conditions for positivity. 

1 Introduction 

The set of positive maps acting on a finite dimensional Hilbert space is a long- 
standing subject of mathematical interest. In spite of many efforts (sec [f 5] 
and references therein) , the structure of this set in spaces of arbitrary dimension 
is still not well understood. Of particular interest are positive maps, which are 
not completely positive [6-8]. The theorem of Jamiolkowski implies [2] that 
any such map can be represented by an operator, acting on a bi-partite Hilbert 
space, which is not positive, but is block-positive. 

Non completely positive maps recently attracted a considerable attention 
of the physics community [9-11]. Positive maps have mainly been studied in 
view of their possible application to characterize quantum entanglement [12] 
and in connection to entanglement witnesses [13-16]. An entanglement witness 
is a Hermitian operator W such that Tr (Wcr) ^ for any separable state a, 
while negativity of Tr (Wp) implies that the state p is entangled. Note that a 
Hermitian operator W may be considered as an observable, so the expectation 
value Tr (Wp) can be measured in an experiment [17]. From a mathematical 
perspective any entanglement witness is a block positive operator which is not 
positive. 

In the present paper we aim to clarify the relation between positive maps and 
positive polynomials. Definitions and basic information can be found in Section 
[2J In Section [31 we explore the link between positive maps and positive poly- 
nomials and we address problems related to early contributions on the subject. 
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In particular, we analyze implications of the work of Jamiolkowski [18, 19] and 
show why the results of these papers do not allow one to formulate a conclusive 
test for positivity of a given map. 

On the other hand, in some particular cases such results can be obtained. In 
Sections 0] and El we investigate two families of maps and working with the cor- 
responding polynomials we find explicit conditions for positivity. Furthermore, 
we demonstrate how positive maps relate to the existence of positive polynomi- 
als which are not sums of squares and we formulate an open problem concerning 
entanglement witnesses in 2 x m dimensional spaces. 

2 Block positivity - motivation and definitions 

Let H\, Hi be finite dimensional spaces over (D, both equipped with Hermitian 
inner products (dim7Yi = N%, dim7Y 2 = A 2 ). Let C {Hi) denote the algebra of 
linear operators on Hi. We denote with C(Hi) + the set of positive elements 
of C{H\). A linear map $ : C{Hi) — > C(Ti.2) is called positive if and only 
if it maps elements of C{H\) + to elements of £(7^2 ) + ■ It is well known [2] 
that the set of positive maps is isomorphic to the set of block positive operators 
(block positive over (D). Therefore, instead of asking whether a given map is 
positive, in this work we will be concerned with the equivalent question whether 
the corresponding operator is block positive, so it can serve as an entanglement 
witness. 

A Hermitian operator A on H — Hi ® H2 is called block positive over (D if 
it satisfies the following condition, 

(u®v\ A(u®v)) > y U £H 1 ,veH 2 - (!) 

Note that condition ([I]) is not invariant with respect to global unitary transfor- 
mations on H, so this definition depends on the particular form of the decom- 
position of H. 

It will also be useful to introduce the concept of block positivity for real linear 
spaces. Let A and Y be finite dimensional vector spaces over 1R (dim A = Mi, 
dim Y = 1/2)- Let A be a linear operator on A (g) Y , In analogy to ([T]), we say 
that A is block positive over 1R if it satisfies 

{x®y) ■ A{x®y)^{) \/ xeX , y eY- (2) 

Condition (T2J) does not imply symmetry of A, but we may always assume that 
A is symmetric because the antisymmetric part of A in (J5|) vanishes. Thence 
(X (g)Y) 3 (wi , W2) 1— > w\ ■ A (11)2) € R is a symmetric bilinear form on X ® Y. 
In index notation, condition ^ reads 

A ab , cd x a y b x c y d > ^ {x a } ^ m fj iCn , (3) 

where x a and y b are the coordinates of x and y with respect to the orthonormal 
bases {e,}^, {fj}^ of A, Y (resp.) which we are using. 

Obviously, §3§ is a positivity condition for a real multivariate polynomial of 
degree 4. If the polynomial A a b tCt ix a y h x c y d is a sum of squares (SOS) of some 
other polynomials Pi, then we must have 

A ab , cd x a y b x c y d = ]T/f = ^(i?>V) 2 , (4) 
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where the real coefficients B l ab (a — 1, . . . , Mi, b = 1, . . . , M%) are arbitrary and 
the range of the index i is finite. 

Indeed, the polynomials Pi must be homogeneous and of degree 2. They 
cannot have terms of the form x a x b , neither of the form y a y b , since there are 
no terms (x a x b ) nor (y a y b ) in the sum A a b, c dX a y b x c y d . Thus we conclude 
that if A a b,cdX a y b x c y d = Yli^i f° r some polynomials Pi, then P t = B ab x a y b . 
But (fj| looks just like a quadratic form on X ® Y, written in the product basis 
{ei ® fj}^{ It is tempting to say that ((4]) implies positive semidefinitness of 
A, but this is not true. Nevertheless, a similar result can be proved if we assume 
that A is symmetric with respect to partial transpose, A T := (1 (g) T) A = A, 
where T denotes the transposition. Putting this in a different way, A should 
satisfy 

(xi ® yi) ■ A(x 2 <%> y 2 ) = (xi ® y 2 ) ■ A(x 2 <8 yi) Vximex^imzy- (5) 

For any operator A being a SOS and expressed by eq. (j4|), we may define the 
following operator A, 




It is easy to see that (x ® y) ■ A (x <8> y) — (x <E) y) ■ A (x <E> y) for all a; G X, 
y £ y. In Appendix El we show that this property together with ([5]) and (J6j) 
imply A = A. But ^4 is of the special form ([5]), which we did not assume about 

A. More precisely, A is proportional to a sum of a semipositive definite operator 
B with matrix elements J^i ^ab^cd an< ^ ^s partial transposition B T with matrix 
elements J2 i B l ad B l cb . We conclude that A a t iC dX a y b x c y d = ^ implies 

A= l -{B + B T ), B^O. (7) 

for the operators A with the property (|5|). A Hermitian operator A is called 
decomposable [1,3] iff A = C + £> T , where C, 13 ^ 0. When ([5|) holds, one can 
easily prove that ([7|) is equivalent to decomposability of A . Thus we arrive at 
the following conclusion, 

Proposition 1. Let X,Y be finite dimensional linear spaces over 1R. Let W 
be the set of block positive, indecomposable operators on X <E>Y which are 
symmetric with respect to transposition and partial transposition. Denote 
with V the set of positive real polynomials of the form A a b^ c dX a y b x c y d which 
are not SOS. There is a linear isomorphism between W and V . 

Proof. The isomorphism in question is II : W 9 A i— > A a b^ c dX a y b x c y d E V . Wc 
still need to show that n is one-to-one. To this end, we assume the equality 
J2a,b,c,d A o.b,cdX a y b x c y d = Y, a ,b,c,d B ab,cdX a y b x c y d for some operators A, B £ 
W. Choose some a, c e {1,2,..., Mi}, b, d S {1,2,..., M 2 }. Considering the 
coefficients at x a y b x c y d in the two polynomials, we get A a b^ c( i+A a( i,cb — B a b iC d + 
B a d,cb- Thanks to the partial transpose symmetry of A and B, we get A — 

B. This tells us that n is injective. On the other hand, every polynomial of 
the form A a }, iC( ix a y b x c y d is an image by II of the partial transpose symmetric 
operator h (A + A T ). The operator | (A + A T ) must be an element of W for 
A a b,cdX a y x c y d to be an element of V (cf. the discussion above). Thus we 
conclude that n is surjective. □ 
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It was demonstrated by Choi [6] and St0rmer [20] that there exist positive 
maps which are not decomposable. The example by Choi can be easily used 
to show that there exist, by Proposition [TJ positive polynomials of the form 
A a b,cdX a y h x c y d which are not SOS [6]. Proposition [T] gives a general motivation 
to investigate block positive operators over R on account of their connection to 
sums of squares. It may also be expedient to study the real case in order to 
develop intuitions about block positivity over <C. It should, however, be kept 
in mind that {1} and (T5J) are not the same. Despite an apparent similarity, 
the block positivity over (D should not be perceived as a simple generalization 
of the block positivity over R. In general both definitions of block positivity 
do not coincide, what can be demonstrated by the following example of a real 
symmetric matrix, 



A = 



This matrix represents an operator on <D 2 ® CD 2 written in the standard product 
basis, {|00), |01), |10), |11)}. It is easy to show that A satisfies inequality (f5]), 
but it does not satisfy condition fl} . Hence the matrix A in (0 is block positive 
over R, but is not block positive over (D. Moreover, when considering operators 
with unit trace, one can easily show that the set of such block positive operators 
over (D is compact whereas block positivity over R does not imply compactness. 
In spite of this basic difference between the two notions of block positivity, 
there exist families of matrices for which conditions and ^ turn out to be 
equivalent - see Section |H 

3 Block positivity and quantifier elimination 

Although the block positivity condition {T]) is simple to understand, it does not 
seem easy to check. The early papers by Jamiolkowski [18,19] suggest that the 
problem can be solved effectively. Even though this conclusion is in some sense 
true, we show a weak point of the argument presented in these papers. 

For convenience of the reader, let us recall the details of the reasoning pre- 
sented in [18]. First, we write condition JT]) in index notation, 

A al3 , 1 su a v l3 u 1 v s ^ V { u * } ^i i){u <s}W2 iCC . (9) 

Next, we introduce blocks, 

{A^) ai := A aPnS v v 5 , (A^) 0S ~ A a ^ lS u a u\ (10) 

We can interpret them simply as matrices or as operators on 7ii and H.2, re- 
spectively. Block positivity condition ^ can be rewritten as 

A { u ] > V ueWl , (11) 

where refers to semipositive dcfmitcness. We shall concentrate on the right 
hand side of (fTT|) . Scmipositivity of A ( u ' is equivalent to the following set of 
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inequalities, 

ueWiVj=i...iV3i (12) 

l^ii<i 2 <— <ii^AT 2 

J is the minor of A> a involving the columns and the rows 

with the numbers ii,...,ii. It follows from the discussion in [18] that the 
functions Wi are homogeneous real polynomials of an even degree in the variables 
{Re(u a )}^ 1 , {Im {u 1 )} 1 ^^. Thus (fT2|) is a set of positivity conditions for real 
homogeneous polynomials of an even degree. If we could solve these conditions 
explicitly, we would answer the question whether a given matrix is block positive. 

That was the idea presented in [18] by Jamiolkowski, who suggested con- 
sidering J2i u i 2 ,...,i n C ixi2-i n X % i X% . . .X%> {Ci li2 ... in S R) as a polynomial in 
the variable X n with coefficients in R[Xl, . . . ,X n _i]. He obtained positivity 
conditions for such a polynomial in a disjunctive normal form, 

V {ii „.„,„_ i} ceV /\D*-(a;i,... > ar n _ 1 )>0 J (13) 

i 3 

where D l - G R \X\, . . . , X n -x] Because the same procedure could be applied 
to any of the Dys (considered as elements of R [Xi, . . . , X n - 2 ] [-X" n _i]), it was 
claimed that the number of variables in (fT3"|) can be iteratively reduced so as to 
yield quantifier free formulas. The problem with this argument is that eq. Q13p 
does not turn out to be equivalent to 

V A ! • •• D i ' ' ■ > °' ( 14 ) 

i 3 

so one cannot use the procedure iteratively. 

To the best of our knowledge, no simple method is known to check positivity 
of a general multivariate polynomial. It is in principle possible to eliminate quan- 
tifiers [21] from formulas like V{ xll ..., Xn _ 1 } C H J2i u i 2 ,...,i n Ci^-in^i ■ ■ ■ x n ^ °i 
but the outcome involves zeros of univariate polynomials of an arbitrary high 
degree, which cannot in general be expressed in terms of the coefficients of the 
polynomials. The known quantifier elimination procedures are laborious and 
should not be expected to provide a constructive solution to the problem. Thus 
we have to conclude this section by repeating the accepted statement that the 
question of explicit conditions for block positivity remains open. 



A three— parameter family of block positive 
matrices 



Fortunately, there exist some particular cases for which positivity conditions 
(TT2")) turn out to be useful in checking block positivity. Let a, b, c € (D. Consider 
the following family of matrices, 



F = 



Foofio -Fbo.oi -Poo, 10 -^00,11 

-Poi.oo -Pbi.oi Fqi,w -Pbi,n 

-Pio.oo -Fio.oi -Pio.io -^10,11 

Fix, oo F\im Fnio Fnii 



2 

a 





a 
i 

2 
h 



(15) 
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which represent operators on Hi ® H-2 — C 2 <8> (D 2 . 

We are going to test condition |T]) using the method suggested in the previous 
section. The blocks (ITU1) with respect to the subsystem described by TL2 are 



Fi 2) {a Ac) 



«l| 2 + |«2| 2 



a \u\\ 2 + c |u 2 | 2 + buiu 2 



a |mi| 2 + c |u 2 | 2 + buiu 2 



Ul\ 2 + \U2\ 2 



(16) 



Obviously, (a, 6, c) is semipositive definite for all u £ <C 2 if and only if 
dctisi 2) (<z,&,c) ^ 0V ueC 2. That is, 

2 



ui 



"2 



a \ui\ 2 + c \u 2 \ 2 + bu\U2 



^ V Uli „ 26 c- 



(17) 



Keeping |ui| and \uq\ constant, we can maximize the absolute value of the term 
a |ui| 2 + c \u2\ +bu\U2 by choosing the phases of u\ and U2 such that the phase 
of bu\U2 is the same as the phase of a|ui| 2 + c|u2| 2 - So done, we see that 
condition (fT7j) is equivalent to 



Q(^ 2 +2/ 2 )) -(\ax 2 +cy 2 \ + \b\xy) 2 ^0 V a , 



J/GIR+- 



(18) 



In inequality (fT5|) , we substituted x for and y for |xt2 1 - It is now easy to see 
that (fT5|) is the same as 



- (x 2 + y 2 ) - \ax 2 +cy 2 \- \b\ xy ^ ^ x , y eTR- 



(19) 



We extended the domain of x, y in (fT!))) to R, which is permissible because \b\ xy 
does not increase if we change the sign of x or y from plus to minus. Substituting 
x — ► r cos ?/ — > r sin ^ in (fT§|) we obtain 



1 - I a + 7 costal - \b\ sin 93 ^ V^m, 



(20) 



where a := a + c, 7 := a — a Condition (|20|) can be easily solved in the two 
following situations: 

a) Re (0-7) = <S=^> \a\ = |c| 

b) Re (07) = ± |a| I7I a = rc, r e 1 
In the case a), condition (|20p simplifies to 



a\ 2 + |7| 2 cos 2 (yO - \b\ sinip ^ V^r. 



(21) 



We observe that \a\ 2 + \j\ 2 ^ 1 must hold in order that (|2ip be true. Keeping 
this in mind, we can rewrite (|2ip as 



A - A + 1 - 



b 


2 






(i«r ^ 




1 







^0 V 



Ae[|a|Vl"l 2 +l7l 2 l 



(22) 
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where we substituted \J\a\ 2 + | 7 | 2 cos 2 ip — > A. As a positivity condition for 
a quadratic function, (f22l) can be easily solved explicitly. Together with the 



condition on \a\ 2 + | 7 | 2 , we obtain 



\a\ 2 + | 7 | 2 < 1 A \a\ + H 2 ^ 1 A |2 |6| 2 |a| < | 7 | 2 V 2 |&| 2 ^M 2 + | 7 | 2 ^ | 7 | S 

which is an equivalent form of (|22[) . In the case b), it is even simpler to get 
the conditions on a, 7 and b equivalent to (f2T)]) . We have \a + 7 cos<y9| ^ \a\ + 
| 7 | |cos . Either for p> or for tp — > n—(p, we obtain \a + 7 cos (p\ — |a|+| 7 | |cos <p\ 
and sin ip is not changed by the substitution <p — > 7r — (/p. Hence we can rewrite 
(f20f as 

1 - \a\ - | 7 | |cos</?| - |6| sinc^ > V^r. (24) 

This is equivalent to (1 — |a| — | 7 | cosp — \b\ sin<p) ^ OV^gn, which is easy to 
solve explicitly in terms of a, 7 and b. We get 



1-M-Vl7| 2 + H 2 ^0. (25) 

In the case of general a, b and c, condition (|20|l is equivalent to the following 
system of four inequalities, 

i) |a| + | 7 | ^ 1, (26) 

ii) M 2 -H 2 < |Re(a 7 )|, (27) 



m 



1-I«| 2 -M 2 



> 2|Re(a 7 )|, (28) 

iv) (| 7 | 2 + |6| 2 ) Z cos 4 V-4(| 7 | 2 + |6| 2 )Re(a 7 )cos 3 V'+ (29) 

+ (4Re (a 7 ) 2 - 2 (l - \a\ 2 - |6| 2 ) (| 7 | 2 + |6| 2 ) - 4 | 7 | 2 ) cos 2 V + 

-4 (^3- |a| 2 - |6| 2 ) Re(a 7 )cosV'+ (l - |a| 2 - |6| 2 ) 2 -4|a| 2 ^ 0, 

which have to be satisfied for all real ip = 2(p. The expression on the left hand 
side of condition iv) is a polynomial P(cosip) of degree four in the variable 
cos ip. This condition means that P is nonnegative in the interval [—1, 1]. Given 
particular values of a, b and c, nonnegativity of P in [—1, 1] can be easily checked 
using the Sturm sequences [22] . It is also possible to produce general conditions 
on a, 6, c in this way, but the resulting formulas would be too complicated to 
reproduce them here and not suitable for further analysis. 

An analogous problem of block positivity over 1R can also be solved for the 
family of matrices (fTS"]) . Most of the work has already been done above. We 
only need to observe that the passage from (flTf to (fTHj) is possible also when 
o, b, c and u\, ui are real numbers. This is true because the maximal value of 
a \ui\ 2 + c |ii2 1 2 + bu\Ui for fixed \u,2\ is \au\ + cu 2 ,] + \b\ \ui\ \v,2\- Thus 
the condition (TTH| turns out to be equivalent to block positivity over 1R of the 
matrices of the form (fTS"]) with a,b,c £ 1R. Later analysis follows as in the case 
b) discussed above. In this way we arrive at two important conclusions. Firstly, 
symmetric matrices of the form (|15[) are block positive over R if and only if 
they are block positive over (D. Secondly, the block positivity condition takes 
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the form (|25[) with a = a + c, 7 = a — c. On the other hand, positivity conditions 
for the family of matrices (fT5|) are easily obtained, 

^-^-^-^ + H 2 N%0Ai-H 2 -|5| 2 -|c|%0. (30) 
We can compare them with the block positivity condition (125p in a picture. 




Figure 1: The grey set of positive semidcfinitc matrices defined by cq. (|30|) is 
contained inside the set of block positive matrices determined by ([23)) . In this 
case the block positivity over (D is equivalent to the block positivity over R. It 
is assumed here that a, b, c £ R, but formulas (fzTl)) and |23|) apply also for a, b, c 
complex, provided that a = rc with r £ R. 



It is clear from Figure [T] that the conditions (f3"0|) and ([23)1 are not equivalent, 
and the set of positive matrices of the family (|15|) forms a proper subset of the 
set of block positive matrices. 

A similar investigation can be performed for a related family of matrices, 



E(s,p,q,r) 



s 
i 

2 
P 







p 

_L 

2 

q 



(31) 



with real parameters s,p, q and r. The block positivity conditions for E (s,p, q, r) 
can be obtained using the methods presented in this section. In particular, tak- 
ing E (a, |, c, |) with a, b, c real, we get a symmetrization of the family (fT5)) . 



F' (a,6,c) := 



F (a, 6, c) + F (a, 6, c) T 



2 

a 



6 
2 



(32) 



Deriving conditions for positivity and block positivity of the matrices F'(a, b, c), 
it turns out that in this case both properties do coincide, unlike in the example 
discussed above. In the light of Proposition [TJ this fact can be understood as a 
consequence of the following theorem [23] , 
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Theorem 2 (Calderon). Let m G N, x 1 , x 2 € R, {V , {^a&cd}a,6, c ,d C R 
Adbcd^xV > 0V x i, x 2 ){j/i}r=i ==* A^xVzV = ^ (B^V) 2 ■ (33) 

TTiai is, anj/ positive biquadratic form in 2 x m variables is the sum of squares 
of quadratic forms □ 

According to Proposition [T] and the Calderon's result, any operator A on 
X®Y = R (8R m which is symmetric with respect to partial transpose and block 
positive over R, is decomposable as well. More than that, we know from the dis- 
cussion preceding Proposition Q] that A can be written in the form (B + B T ) /2 
with B ^ 0. In the case of A = F' (a, b, c), B must be of the form (f5T|) with 
s = a, q = b and p + r = b. As can be checked by direct computation, the char- 
acteristic polynomials of E (s,p, q,r) and E (s,r,q,p) are the same. It follows 
that E <^> F T ^ 0, which in turn leads us to the conclusion that the matrix 
F' = | (F + F T ) = \ (E + E T ) is block positive if and only if it is positive. 

To explain our observation about F', we could also have used the St0rmer- 
Woronowicz theorem [1,3], which implies that an operator A on R 2 ® R 2 (or 
on R 2 (g) R 3 J is block positive if and only if it is decomposable. This suggests a 
possible connection between the Calderon's and the St0rmer-Woronowicz the- 
orems. On the other hand, the first theorem holds for all R 2 (g> R m (m G N) 
whereas the latter works for m < 3 only. 

The theorem of Calderon allows us to find some further implications for the 
subject of positive maps. 

Proposition 3. Let m G N. Either all block positive operators on (D 2 (g> <C m 
with real matrices are decomposable or there exists an operator A on (D 2 ® 
(D™ 1 with real matrix elements A a b yC d such that A a b^ c dX a y h x c y d is the sum of 
squares of bilinear forms, but A is not decomposable. 

Proof. Let A be an operator on (D 2 ® (D™ 1 with real matrix elements. If A is block 
positive on (D 2 ® CI m , it must be block positive on R 2 ® R m . From Calderon's 
theorem it follows that A a b^ c dX a y b x c y d is the sum of squares of bilinear forms. If 
this implies decomposability of A, any block positive operator on (D 2 (C m with 
real matrix elements is decomposable. If not, there exists an indecomposable 
operator A such that A a b^ c dX a y b x c y d is SOS. □ 

Both the mutually exclusive possibilities in Proposition[3]are interesting and 
it will be good to know which of them is true for which m (of course, the answer 
is known for m = 1, 2, 3 - every positive map is decomposable). We hope that 
stronger results of similar kind can also be obtained and they should give better 
insights into the structure of positive and indecomposable maps. 



5 Block positivity of 4 x 4 matrices over R 

We want to illustrate the abstract discussion presented in Section [3] with a 
concrete example. To that aim, following [25], we derive sufficient and neces- 
sary conditions for an arbitrary operator A on R 2 ® R 2 to be block positive 
Let the matrix elements of A be A a f, tC d (a, b, c, d G {1, 2}). The blocks with re- 
spect to the first subsystem have the matrix elements [A^ J = A a ± >c i (y 1 ) + 
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(Aai,c2 + A a2 .ci)y 1 y 2 + A a2 , C 2 (y 2 ) 2 - Positivity of Ay is equivalent to the re- 
quirements that Tr Ay ^ and det Ay 1 ^ 0. Nonnegativity of the trace of 
Ay 1 ^ for all y = (j/ 1 ,?/ 2 ) G R 2 means that 

2 2 2 

X!^mi (?/) + X! + A ]2 , 3 i)y 1 y 2 + ^Afe2,fc2 (y 2 ) ^0 V y i, a 2 eR . 
i=i j=i fe=i 

(34) 

Obviously, (I34j) is a positivity condition for a quadratic form on R 2 and we can 
write it explicitly as 

2 2 2 if 2 

ij — l i—1 k—1 

The expression for the determinant of A y reads 

det A y 1] = c 4 x 4 + c 3 x 3 z + c 2 x 2 z 2 + axz 3 + c z 4 , (36) 
where we substituted x for y , z for y 2 and we introduced 
co = ^1242^21,21 — ^4.ii,2i^4.2i,ii) (37) 

Cl = ^4.22,22 (^4.12,11 + ^4.11,12) + ^-12,12 (^-22,21 + ^-21,22) + (38) 

— ^■22,12 (^-12,21 + ^-11,22) — j4ll,21 (^22,11 + -421,12) , 

C2 = -All, 11^22,22 + ^.21,21^-12,12 + (A-11,12 + 4.12, n) (A.21,22 + A 22 ,2l) + (39) 

— ^.11,21^.21, 11 — Ai2, 22^.22, 12 — (A-11,22 + -4i2,2l) (-4.21,12 + ^22, ll) j 
C.3 = -4-11,11 (^21,22 + ^22,2l) + (^-11,12 + Al2,u) + (40) 

— -4.11,21 (A-21,12 + -4-22, ll) — A-12.YI (-4n, 22 + -4l2,2l) , 

C 4 = An, 11-421, 21 — -4.11,21-4.21, 11- (41) 

The Cj's are homogeneous polynomials in the matrix elements 4. a b. c d. It is easy 
to see that non-negativity of (f3"B"|) for all x, z G 1R is equivalent to 

c 4 :r 4 + c 3 .T 3 + c 2 a; 2 + cix + c > V KeR . (42) 

Thus we showed that in the case of a symmetric matrix A of order 4 condition 
© is equivalent to ([33)1 plus (|4"2"]> . The inequalities ([55)) are explicit conditions 
on the matrix elements A a b iC d, but in ([4"2")l we need some additional work to 
dispose of the quantifier v^r. There is no single method of doing it, but the 
one which seems most economical to us is by using the following theorem [24] , 

Theorem 4 (Sturm). Let f = fo be a real univariate polynomial with no mul- 
tiple roots in R. Let f\ be the first derivative of f . Define 

/„+! := rem (/„_i, /„) , (43) 

where rem(/i,g) is the remainder obtained when dividing h by g. Define N (r) 
as the number of sign changes in the sequence 

fo (r) , fi M , ~h (r) , -h (r) , h (r) , h M , ~h (r) , . • • (44) 
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with zeros skipped. Assume a, (3 € R, a < (3, fo (a) ^ and fo (/?). TTie number 
of zeros of fo in the interval (a, f3) equals N (a) — N {(3), 

N(a)~N (/?) = # {r e (a; /?) \f (r) = 0} . (45) 

□ 

The sequence of functions (|43[) is the same as in the Euclid's algorithm 
applied to / and /'. When the signs are changed as in (|4"4"|) . the sequence is 
called the Sturm sequence of /. We know that {/n} Tl=0 i must terminate at 
some f m £ R\ {0}, which is the greatest common divisor of / and /'. If we go 
to the limits a — — oo, (3 = +oo in Theorem |4j we easily obtain the number of 
real roots of /, 

Corollary 5. Let f — fo be a real univariate polynomial with no multiple roots 
in R and f\ - its first derivative. Let f n (n = 2,3...) be defined like in |^3| ) 
and assume 

fn (r) = a ntkn r k " + a n ,k n -ir 1 ' n ~ 1 + . . . + a ,„, (46) 

where k n 0, a n ,fc n ^ V n . Denote with N (+oo) the number of sign changes 
in the sequence 

O>0M) J a l,fcl i ~ a 2,fc 2 ' ~ a 3M 7 a 4,fe 4 7 ■ • ■ ' ± a m,fc m (47) 

and with N (—00) the number of sign changes in 

I— J a o,fe )l — J a l,fei)l — J a 2 ,fc 2 , v J a3,fc 3 7l~J a i ^ kl , . . . ,±a m .k m . 

(48) 

T/ie number of real zeros of f equals N (—00) — N (+00), 

N (-00) - N (+00) = # {x e R|/ (x) = 0} . (49) 

□ 

Let us take for / the polynomial C42; 4 + C3X 3 + c^x 1 + C1X + C0 which appears 
in (|42p. We shall now assume that it has no multiple roots in R. Then we can 
use Corollary to check positivity of /. 

The sequence {fn} n= o 1 consists of at most five polynomials, 

f = fo = C4X 4 + C3.T 3 + c 2 x 2 + c\x + Co, (50) 
fx = 4c 4 x 3 + 3c 3 a; 2 + 2c 2 x + c\ , 
h= CL2.1X 1 + a 2 .\x + a 2 ,o, 

h= a3,l ;E + a 3,0, 

fi = a4,0- 

If we make an additional normality assumption, which says that the degrees 
of fo , . . . , fi drop one by one in the successive lines of (|50p , it is easy to write 
down positivity conditions for /, 

c 4 > A (o 2 ,2 > V 013,1 > 0) A a 4:0 > 0. (51) 

The expressions for 02,2, 03,1 and 0,4.0 can also be easily obtained in the present 
situation. We get 



1 32c 4 » 

02,2 = -r- — 01, a 3 1 = — 0-0-2, a.4,0 = -771 2 a 3> l 52 J 

16c 4 erf 64C40J 
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where 

a i := 8C2C4 — 3c§, (53) 

(72 := 3ciCg — I4C1C2C3C4 — c| (c| — 6C0C4) + 2c4 (2c 3 , + 9c 2 C4 — 8C0C2C4) , (54) 
03 := 2ci (2ci 2 - 9c c 2 ) c 3 3 + 2cic 3 c 4 (-9c 1 2 c 2 + 40c c 2 2 + 96c 2 c 4 ) + (55) 
27c 2 c 3 4 + c 3 2 (-ci 2 c 2 2 + 4c c 2 3 + 6c ci 2 c 4 - 144c 2 C2C 4 ) + 
c 4 (4ci 2 c 2 3 + 27ci 4 c 4 + 128co 2 c 2 2 c 4 - 256c 3 c 4 2 - 16c c 2 (c 2 3 + 9ci 2 c 4 )) . 

According to ()52l) . the normality assumption is equivalent to C4 ^ A a± 7^ 
A (72 7^ A 0-3 7^ 0. If these conditions hold, we can rewrite (f5Tj) as 

c 4 > A ((7i > V cr 2 > 0) A (7 3 < 0. (56) 

This is an explicit condition for / to be positive. Of course, it was obtained 
under the assumption that {/n}^ =0 are normal. Nevertheless, we can use (1561) as 
a starting point for an all-purpose nonnegativity test for polynomials of degree 
less or equal four. Indeed, suppose that the sequence {f n }n=o 1S n °t normal. 
That is, at least one of the numbers C4, a\, (72, (73 happens to be zero. If C4 = 0, 
the nonnegativity question becomes trivial. We get that C3X 3 + c 2 x 2 + C\x + cq 
is nonnegative if and only if 

c 3 = A c 2 > A c 2 - 4c 2 c < 0. (57) 

The case in which C4 7^ but 01 02(73 = can be analyzed using a little more 
sophisticated techniques (see Appendix |B|). All in all, we arrive at the following 
nonnegativity conditions for /, 

{c 4 > A ((7i ^ V (7 2 ^ 0) A cr 3 <0} V 

V {c 4 = A c 3 = A c 2 A c\ - 4c 2 c < 0} , (58) 

where (73<0 means 3^>oV^'<^ ((73 (04, C3, C2, c\, cq + £') < 0). We can write <73<0 
explicitly as 

(73 < V (03 = A (ki < V (ki = A n 2 sC 0))) , (59) 

where 

Ki = 4c 3 c — I8C3C2C1 + 80c4C3C 2 ci + 6c±c 3 c 1 — 16c4C 2 (c 2 + 9040]^) , (60) 
k 2 = 27c 4 - 144c 4 c^c 2 + 128c^c 2 + 192e|c 3 Ci. (61) 

Obviously, conditions |58|) and (|35|) together with the definitions (|59|). (|60|) . 
(|6T|) . (pj) . (|54l). ([55]) . (137]) . (l38|) . ([39]) . 140]) and gl]) provide us with a method 
to test block positivity over 1R of 4 x 4 matrices. We see that lengthy calculations 
are involved, even though the studied example is the simplest possible one. It 
is also clear that the iterative procedure proposed in [18] could not work with 
conditions like (|56|) . let alone (fS"8"l) . 

6 Conclusions 

We have re-examined the method [18] of establishing positivity of a map with 
help of multivariate polynomials and we conclude that in the general case this 
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problem remains open. The same can be said about the equivalent problem of 
checking whether a given operator acting on a composite Hilbert space is block 
positive. Nevertheless, for certain family of operators checking the positivity 
of the associated polynomials allowed us to find concrete criterion for block 
positivity. Such concrete examples are provided in Sections |4] and [5] By giving 
the example ([5]). we touched upon the relation between the block positivity 
conditions over (D and over R. 

We also outlined connections between block positivity, indecomposability 
and the sums of squares (Propositions Q] and [3J Theorem \2fy ■ Proposition [3] 
opens a discussion about the two mutually exclusive possibilities concerning 
indecomposable maps on (D 2 ® (D m (cf. Section |4]). 

Finally, we tried to show that polynomials, which have been thoroughly stud- 
ied by mathematicians and engineers, may deserve more respect of physicists 
working on quantum information or on open quantum systems. In particular, 
the separability problem itself can be formulated as a set of polynomial equal- 
ities [26]. Techniques like the calculation of a Grobner basis of an ideal are 
widely used to solve polynomial equations and they could be of importance in 
physical problems like the separability problem. 
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cocos. 

A Equality between operators A and A 

Consider the symmetric bilinear forms $ : (X (8 Y) 2 3 (w\, W2) 1— > wi ■ A (11)2) £ 
R, # : {X(g>Y) 2 3 (wi,w 2 ) h-> wi ■ A(w 2 ) £ R. We know that $(x(g)j/) = 
$ (x (g> y) for arbitrary x £ X, y £ Y. From |5]) and (J6j> it follows that <£>, $ are 
symmetric with respect to partial transposition, 

$ (xi ® yx,x 2 <8> y 2 ) = $ (xi <g> y2,%2 , (62) 
®(xi ®y\,x2 ® 2/2) = * (xi ® 2/2,2:2 ® yi) ■ (63) 

Let us choose x £ X and consider the following maps, & x : Y 2 3 y 1— > 
$ (xg> j/i,x® y 2 ) G R, $2 : ^ 2 3 2/ >-> $ (a: ® 2/1, x ® y 2 ) G R. From 
and (I63p we know that $3;, $3; are symmetric bilinear forms on Y. As a con- 
sequence of $ (x ® y) = $ (x ® y), $2, (y, y) = $3, (y, y) for arbitrary y £ Y. 
Hence the quadratic forms corresponding to $ and $ are equal. This implies 
$ x = so we get 

$ (x (g) yi,x «) y 2 ) = i(x® j/i,x® j/2) Vxex, yi ,y 2 eY- (64) 

Now we consider the maps $j, 1)1(2 : A 2 9 (xi, x 2 ) 1— > $ (xi ® j/i, #2 ® J/2)) &yi,y 2 '■ 
X 2 3 (xi,x 2 ) 1— > $ (xi ® J/i,X2 <8> j/ 2 ). From the symmetry of $, $ and the 
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properties (|62[) , (|63p . we see that <& yi , y2 , & yi , y2 arc symmetric bilinear forms on 
X. As a consequence of (|64|) . $ yi , y2 (x, x) = & yi , y2 {%,%) for all a; S X. This 
implies & yi , y2 (xi,X2) — Q yi .y 2 (^1,^2) for arbitrary 2:4, ir 2 G X. In this way 
we get $ (xi (g> yi, £ 2 <8> 2/2) = $ (^l <8> 2/1, £ 2 ® 2/ 2 ) V^xaex.^^er, which is the 
same as 

(xi ® yi) • A (a; 2 O 2/2) = (ai <8> 2/4) • A (x 2 <8> 2/2) ^x u x 2 ex^ yuy2 er- (65) 
Of course, implies A = A. 

B Nonnegative polynomials with <7i0"20"3 = and 

C4^0 

Our aim is to figure out all the sign configurations of C4, 01, 02, 03 such that they 
meet the constraints C4 5^ A 040203 = and they correspond to nonnegative 
polynomials / = C4X 4 + c 3 x 3 + c 2 x 2 + ax 1 + cq. We also have to check that 
the remaining sign configurations can never give a nonnegative /. Of course, 
c 4 < implies that / (x) be negative for some x, so we only need to consider 
C4 positive. First we show that 03 > cannot happen for a nonnegative /. 
Suppose 03 > 0. We know that 04 = or er 2 = 0. Let us first consider 04 = 0. 
Because 04 = 8c 2 C4 — 3c| and C4 > 0, we can increase c 2 by £ > and get eri > 
for sure. If c 2 turns out to be zero after this operation, we additionally increase 
Co by £ > 0, which must give us ct 2 ^ because cr 2 = ct 2 — 2C0C4CT4 where <r 2 does 
not depend on cq- The numbers e, £ can be made arbitrarily small, so as not 
to influence the sign of 0-3 . Hence we see that / + ex 2 + £ has a normal Sturm 
sequence and it does not satisfy (l56|) because the 03 corresponding to f + ex 2 + £, 
is positive. But / + ex 2 + £ / implies / ^ 0, so / cannot be nonnegative. We 
conclude that / ^ is impossible for C4 > 0, 03 > and 01 = 0. For C4 > 0, 
03 > and 0" 2 = 0, we only need to increase / by a sufficiently small £ to get to 
the conclusion / ^ 0. Our observations mean that 03 > always implies / ^ 0. 
Let us now consider the polynomials / for which conditions 

c 4 > A (01 = V cr 2 = 0) A 03 < 0. (66) 

are satisfied. If 01 vanishes, we can get 04 > by increasing c 2 (c 2 — * c 2 +e). If 
ct 2 turns out to be afterwards, any change of Cq (cq — > Co+£) will give us cr 2 7^ 
(cf. the discussion above). We can take e and £ arbitrarily small, which allows us 
to avoid changing the sign of 03. After all, we get a polynomial J + ea; 2 +£ which 
has a normal Sturm sequence and it is positive since C4 > A 04 > A 03 > 
for the corresponding C4, 04 and 03. Because e and £ can be arbitrarily small, 
we see that / is a pointwise limit of a sequence of positive polynomials. Hence 
/ is nonnegative. The same conclusion can be drawn for C4 > 0, 04 7^ 0, 0" 2 = 
and 03 < 0, so we should add (p^]) to our set of non-negativity conditions. We 
can write ((MJ) and (|56|) as a single condition, 

c 4 > A (04 ^ V o- 2 ^ 0) A cr 3 < . (67) 

The only situation which is left to analyze is that of 0-3 = 0. To that end, let 
us write 03 as a polynomial in cq, 

03 = K3C0 + « 2 Co + K4C0 + Ko, (68) 
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where 



k = ~c\c\c\ + \c±c\c\ + iclcl - 18c 4 c 3 c 2 Ci + 27c\c\, (69) 
«i = 4cgc 3 - I8C3C2C1 + 8OC4C3C2C1 + 6C4C3C1 - 16c 4 c 2 (C2 + 9c 4 Ci) , (70) 
k 2 = 27^ - 144c 4 c^c 2 + Y2&c\c\ + 192c|c 3 ci, (71) 

K3 = -2564 ( 72 ) 

Because of the assumption c 4 > 0, we know that 03 is not constant with respect 
to Co- The idea now is to infinitesimally increase cq and see what the outcome 
is. If 03 becomes positive, we conclude that / ^ 0. If it turns out to be negative 
(we denote this with /<0), we go back to the initial values of Cj and ask about 
the signs of o~\, 02. If <J\ > 0, we choose £ > so small that <j\ > holds when 
we increase cq by £. Then (j67|) is true for / + £, so / + £ is nonnegative. Since 
the £ in / + £ can be made arbitrarily small, we get / ^ 0. If o\ = 0, we increase 
c o by £ to get (73 < and then we increase C2 by e so as to get o~\ > and not 
to violate (73 > 0. After that ((671) holds for / + ex 2 + £ and again we get to the 
conclusion that / 0. Therefore we can add 

c 4 > A <ti > A ct 3 <0. (73) 

to our list of non-negativity conditions for /. Now we only need to analyze 
the case c 4 > A o\ < A (72 ^ A(7 3 <0 to finish our work. If (7 2 > 0, we 
choose £ > so small that the 02 corresponding to / + £ is also positive. Then 
/ + £ > and we get / > 0. The case a 2 — is also simple to analyze. Because 
<7i < 0, increasing cq causes <7 2 = — c 4 co<7i to become positive, so we get 
<7 2 > A (73 < for / + £ and again this leads us to / ^ 0. Thus we can add 

c 4 >0 A (7i < A (7 2 ^ A (7 3 <0. (74) 

to our nonnegativity conditions for /. It is convenient to write (|67p and (|73|) in 
a single formula, 

c 4 > A ((7i ^ V cr 2 ^ 0) A cr 3 <0. (75) 

Because of the particular form (|68| of (73, we explicitly write the condition <7 3 <0 
as 

cr 3 < V (ct 3 = A (ki < V (ki = A k 2 < 0))J , (76) 
since condition c 4 > implies that K3 < 0. 
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